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Abstract. We present the results of a number of high resolution radio observations of the AGN 1156-1-295. These 
^ ' include multi-epoch and multi-frequency VLBI, VSOP, MERLIN and VLA observations made over a period of 

' 50 months. The 5 GHz MERLIN images trace a straight jet extending to ~ 2" at P. A. ~ -18°. Extended low 

brightness emission was detected in the MERLIN observation at 1.6 GHz and the VLA observation at 8.5 GHz 
with a bend of ~ 90° at the end of the 2 arcsecond jet. A region of similar diffuse emission is also seen about 2 
' - ' , arcseconds south of the radio core. The VLBI images of the blazar reveal a core-jet structure with an oscillating 

' jet on a milli- arcsecond (mas) scale which aligns with the arcsecond jet at a distance of several tens of milli- 

I arcseconds from the core. This probably indicates that the orientation of the jet structure is close to the line of 

■~>^ ■ sight, with the northern jet being relativistically beamed toward us. In this scenario the diffuse emission to the 

i-C ' north and south is not beamed and appears symmetrical. For the northern jet at the mas scale, proper motions of 

13.7 ± 3.5, 10.6 ± 2.8, and 11.8 ± 2.8 c are measured in three distinct components of the jet (go = 0.5, Ho — 65 km 
Q ■ s~^Mpc~^ are used through out this paper). Highly polarised emission is detected on VLBI scales in the region 

5-H ' in which the jet bends sharply to the north-west. The spectral index distribution of the source shows that the 

, strongest compact component has a flat spectrum, and the extended jet has a steep spectrum. A helical trajectory 

' along the surface of a cone was proposed based on the conservation laws for kinetic energy and momentum to 

explain the observed phenomena, which is in a good agreement with the observed results on scales of 1 mas to 1 



arcsec. 
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1. Introduction 

The AGN 1156-1-295 is extremely variable over a 
broad range of the electromagnetic spectrum, from ra- 
dio waves to 7-rays. At a redshift of z = 0.729 
HVeron-C etty fc Veron 19981 ), it has been classified as 
both a Highly Polarised Quasar (HPQ) and an Optically 
Violent Variable (OW ) source IjWills et al. 19831 1992, 
IGlassgold et al. 1983D . In its active phase, the source 
shows fluctuations at optical wavelengths with an am- 
plitude of ~ 5 — 7% on a time scale of 0.5 hour 
II Wills et al. 198311 . There are also large and rapid changes 
of the optical linear polarization, with the percentage 
of polarised flux density changing between 1 — 29% 
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IjWills et al. 1992ll . It is one of the class of radio sources 
detected by EGRET (Energetic Gamma Ray Experiment 
Telescope) on the Compton Gamma Ray Observatory. 
While quiescent 7-ray emission remains undetected in 
the source, three strong flares (^ 3.9 x 10~^° Jy) were 
detected by EGRET at energies > 100 MeV during 
the period from 1992 to 1996 ([Thompson et al. 1995[ 
IMukherjee et al. 1997| IHartman et al. 1999|l . 

Radio images of 1156+295 show a typical 'core - 
jet - lobc(s)' morphology. On the arcsecond scale, the 
VLA 1.4 GHz image shows that the source has a sym- 
metrical structure elongated in the north-south direction 
dAntonucci fc Ulvestad, 1985| ). Th e MERLIN image at 
1.6 GHz and VLA image at 5 GHz (McH ardy et al. 1990| ) 
show a knotty jet extending to about 2" north of the core 
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Table 1. The epochs, frequencies, and arrays of the observations described in this paper 



Obs. 


Epoch 


Band 


Bandwidth 


On-source 


Array 


Telescopes" 


Correlator 


No. 




(GHz) 


(MHz) 


(hrs) 








1 


1996.43 


5.0 


64 


0.06 


VLBA 


AU 10 


VLBA" 


2 


1997.14 


5.0 


28 


10 


EVN + 


Ef Sh Cm Jb Mc On Hh Ur Wb Tr 


MPIfR'^ 






5.0 


14 


10 


MERLIN 




MERLIN 


3 


1997.41 


1.6 


14 


10 


MERLIN 




MERLIN 


4 


1997.42 


1.6 


32 


2.7 


VSOP 


HALCA + VLBA (aU 10) 


VLBA 


5 


1998.12 


5.0 


32 


2.5 


Global 


VLBA Ef Sh Jb Mc Nt On Wb Tr 


VLBA 


6 


1999.01 


15.0 


64 


0.5 


VLBA 


AU 10 


VLBA 


7 


1999.14 


5.0 


28 


10 


EVN + 


Ef Sh Cm Jb Nt On Hh Ur Wb Tr 


MPIfR 






5.0 


14 


10 


MERLIN 




MERLIN 


8 


1999.45 


5.0 


32 


2.5 


Global 


VLBA Ef Jb Mc Nt On Wb 


VLBA 


gd 


2000.15 


1.6 


32 


0.5 


VLBA 


AU 10 


VLBA 


10 


2000.92 


8.5 


50 


0.15 


VLA 


AU 27 


VLA 






22.5 


50 


0.15 


VLA 


AU 27 


VLA 



" Telescope codes: Ef: Effelsberg, Sh: Shanghai, Cm: Cambridge, Jb: JodreU Bank (MK2), Mc: Medicina, Nt: Noto, 

On: Onsala, Hh: Hartebeesthoek, Ur: Urumqi, Wb: WSRT, Tr: Torun 
^ The NRAO VLBA correlator (Socorro, USA) 

The MKIII correlator at MPIfR (Bonn, Germany) 

Polarization observation 



in P. A. = —19°. Diffuse emission is observed both to the 
north and south of the core. 

High resolution VLBI images of the source show 
a 'core - jet' structure with an apparently bent 
jet (p cHard y et al. 1990| [Finer fc Kingham, 19971 
[Hong et al. 1999| . A wide range of superluminal veloci- 
ties of the VLBI jet components have been reported, up 
to 40c (McHardy et al., 1990, 1993). Finer & Kingham 
(1997) reported a slower superluminal velocity in the 
range of 5.4 to 13.5 c, based on 10 epochs of geode- 
tic VLBI observations. Jorstad et al. (2001) reported 
superluminal motion in the range of 11.8 to 18.8 c on 
the basis of 22 GHz VLBA (Very Long Baseline Array) 
observations. Kellermann et al. (1999) and Jorstad et 
al. (2001) find that 7-ray loud radio sources are more 
likely to exhibit extreme superluminal motion in their 
radio jet components. The data also indicate that radio 
flares in 7-loud sources are much stronger than those 
for 7-quiet systems. In highly variable radio sources, 
7-ray outbursts are often associated with radio flares 
UTornikoski fc Lahteenmaki 2000|l . 

In this paper, we present the results of MERLIN ob- 
servations at 1.6 and 5 GHz, VLA observations at 8.5 
and 22.5 GHz, and VLBI observations conducted with 
the EVN (European VLBI Network), VLBA and VSOF 
(VLBI Space Observatory Programme), at frequencies 
of 1.6, 5, and 15 GHz. The images obtained allow us 
to present the morphology of the source from parsec to 
kilo-parsec scales. We discuss how the radio morphol- 
ogy and details of the jet structure can be explained in 
the framework of the standard relativistic beaming model 
iPandford fc Rees 1974,1 . 



2. The observations and data reduction 

The epochs, frequencies, and arrays of various ob- 
servations described in this paper are summarized in 
Table The overall composition of the data sets listed 
in Table ^ consists of observations proposed specif- 
ically for this project (observations #2 fc #7), to- 
gether with additional data made available to us by 
other observers (observation #1 - IFomalont et al. 2000L 
#2 - |Hong et al. l999| #3 - M ERLIN observation of 
a phase calibrator, #4 - [Hirabayashi et al. 1998, #5 - 
m - K^ur vits et al. 20031 ^ 8 - 
#9 & #10 - IHong et al."2003| ). 



Garrington et al. 1999 



Garrington et al. 2001 



2.1. MERLIN observations 

Two full-track observations (EVN + MERLIN) at 5 GHz 
were carried out on 21 February 1997 and 19 February 
1999. Another observation (in which 1156-1-295 was used 
as a primary phase calibrator) was made at 1.6 GHz on 
28 May 1997. The MERLIN array consists of 6 anten- 
nas (Defford, Cambridge, Knockin, Darnhall, MK2, and 
Tabley) at 5 GHz and 7 antennas (the same as above plus 
the Lovell telescope) at 1.6 GHz (jThomasson 1986)l . 



2.2. VLA observations 

The source was observed with the VLA in A-configuration 
at 8.5 and 22.5 GHz as a part of a sample of EGRET- 
detected AGNs in December 2000. All 27 antennas partic- 
ipated in the observation. 



X. Y. Hong et al.: A relativistic helical jet in the 7-ray AGN 1156+295 



3 





Dl 





O 



O 




o 



Fig. 1. MERLIN 5 GHz images of 1156+295 at the epochs 1997.14 (a) and 1999.14 (b), and the 22.5 GHz VLA image 
at the epoch 2000.9 (c). The images were restored with a 80 mas FWHM Gaussian beam. Peak flux densities are 1.6, 
2.2 and 1.6 Jy/beam, the r.m.s. noise values are 0.15, 0.15, and 0.3 mJy/beam, and the lowest contours are 0.5, 0.5 
and 1.0 mJy/beam, respectively. Contour levels increase by a factor of 2. 



2.3. VLBI observations 

Two epochs of full track observations with EVN (+ 
MERLIN), two epochs of Global VLBI (1156+295 again 
was used as a phase calibrator), and three epochs of VLB A 
observations of 1156+295 were carried out from 1996 to 
2000 (see Table 1). The polarization VLBA observations 
(epoch 2000.15) at 1.6 GHz were carried out for 24 hours 
for a sample of EGRET-detected AGNs. 

The EVN data were correlated at the MKIII processor 
at MPIfR (Bonn, Germany), the VLBA and Global VLBI 
data were correlated at the NRAO (Socorro, NM, USA). 

2.4. VSOP observation 

The source 1156+295 was observed by VSOP at 1.6 GHz 
on 5 June 1997 as part of the VSOP in-orbit check-out 
procedure ( |Hirabayashi et al. 1998| ). In that observation, 
ground support to the HALCA satellite was provided by 
the Green Bank tracking station, and the co-observing 
ground-based array was the VLBA. 

2.5. Data reduction 

The MERLIN data were calibrated with the suite of D- 
programs IjThomasson 1986jl . The VLA data were cali- 
brated in AIPS (Astronomical Image Processing System) 
using standard procedures and flux density calibrators 



Table 2. Flux density values adopted for the VLA pri- 
mary calibrators 



Calibrator 


X band (8.5 GHz) 


K band (22.5 GHz) 




(Jy) 


(Jy) 


3C84 


3.22 


1.17 


3C286 


5.18 


2.50 



3C48 and 3C286. Their adopted flux density values are 
listed in Table [3 

The VLBI data were calibrated and corrected for resid- 
ual delay and delay rate using the standard AIPS analysis 
tasks. For the polarization observations (epoch 2000.15), 
the solution for the instrumental polarization (D-terms) 
was based on 16 scans of the calibrator DA193, covering 
a large range of parallactic angles. The absolute polariza- 
tion angle of the calibrator, DA193, was assumed to be 
-58° based on the VLA measurements in D configuration 
at 1.67 GHz on 17 July 2000 (C. Carilli, private communi- 
cation) . This was the measurement closest in time of the 
absolute angle available, 0.4 years apart from the epoch 
of our observation. 

Post-processing including editing, phase and ampli- 
tude self-calibration, and imaging of the data were con- 
ducted in the AIPS and DIFMAP packages (Shepherd et 
al. 1994). The final images were plotted within AIPS and 
the Caltech VLBI packages. 
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The task MODELFIT in the DIFMAP program was 
used to fit models of the source structure. This con- 
sisted of fitting and optimising a small number of elliptical 
Gaussian components to the MERLIN, VLA and VLBI 
visibility data. 

We re-imaged each pair of data-sets with the same m, v- 
ranges, cell size, and restoring beams to produce differen- 
tial images and a spectral index distribution. Each pair of 
data-sets is selected from the same correlator and as close 
as possible in time to extract spectral index data and to 
deduce the offset. 

The data used in reconstructing the spectral index 
distributions were obtained non-simultaneously (except 
the VLA observations at 8.5 and 22.5 GHz). In princi- 
ple, structural variability (proper motion and/or changes 
of size or brightness of the components) will mask the 
true spectral index distribution. In order to minimize 
this masking effect we align images obtained at different 
epochs and frequencies at the position of the core assum- 
ing its opacity shift is small comparing to our angular 
resolution. 

3. Results 

3.1. The arcsecond-scale images 

The two 5 GHz MERLIN images and the 22.5 GHz VLA 
image restored with the same circular 80 mas beam are 
shown in Fig. 1. The peak brightness at 5 GHz has in- 
creased by 37% from 1.6 Jy/beam at the epoch 1997.14 to 
2.2 Jy/beam at the epoch 1999.14. This could be related 
to an outburst in the radio core (see section 3.3). 

The arcsecond-scale morphologies are similar with an 
almost straight jet at a position angle of about —18° with 
perhaps some evidence of a sinusoidal fluctuation on the 
sub-arcsecond scale (Fig. 1). A knot at ~ 0.7" and a 
hotspot 2" from the core are well detected. Within 1" from 
the core, the jet is resolved into several regularly spaced 
knots. When the jet passes the knot, it bends slightly 
and ends with a hotspot around ~ 2" from the core. 
No counter-jet emission is detected with the MERLIN at 
5 GHz and VLA at 22.5 GHz. Only the three main discrete 
components (D, D3 and Dl) were detected at 22.5 GHz 
(Fig. Ic). 

The MERLIN image at 1.6 GHz and the VLA image at 
8.5 GHz are presented in Fig. 2. We restored the two data 
sets with a circular Gaussian beam of 250 mas in diameter 
for comparison. They both show a 2" jet at P. A. ~ —18°. 
Low brightness extended emission was also detected with 
MERLIN at 1.6 GHz and the VLA at 8.5 GHz. This emis- 
sion is seen to bend away eastwards 90°) from the 
main arcsecond scale jet structure. A region of diffuse ra- 
dio emission is seen about 2" south of the core too. 

One explanation for the overall morphology of the 
source is that it is a double-lobed radio source seen al- 
most end-on with the northern jet relativistically beamed 
towards us. Doppler boosting makes the northern jet much 
brighter than its de-boosted southern counterpart. The 



southern jet remains undetected except when it becomes 
sub-relativistic at the end of the jet. It appears as the 
diffuse emission seen to the south of the core. 

Another possible explanation of the brightening of the 
southern jet at distances of 1.4 and 2.2 arcseconds from 
the core is that it bends closer to the line of sight, thus in- 
creasing the Doppler boosting. However, this explanation 
seems to be less likely since it would require the southern 
jet to be closer to the sky plane in its inner parts while 
the inner northern counterpart appears to be close to the 
line of sight. 

3.2. VLBI images 

VLBI images of 1156-1-295 at various frequencies (1.6, 5, 
and 15 GHz) are shown in Fig. 3 in the time sequence from 
(a) to (h). The parameters of the images are summarized 
in Table El 

The eight VLBI images of the source all show an oscil- 
latory jet structure on mas scales. The jet initially points 
almost exactly to the north but then bends to the north- 
east at 3~4 mas from the core. Several tens of milliarcsec- 
onds from the core it finally turns about 90° to the north- 
west, thus aligning with the direction of the arcsecond- 
scale jet. 

The five epochs of VLBI observations at 5 GHz from 
1996 to 1999 allows us to study the structural varia- 
tions of 1156-1-295 in detail. The two global VLBI images 
(Fig. 3 d&g) both with excellent (u, t')-coverage and high 
sensitivity, enable us to detect features that have relatively 
low surface brightness. 

The highest resolution image was made with the 
VLB A at 15 GHz. An additional component has been de- 
tected in the north about 1.5 mas from the core (Fig. 3e). 

The resolution of the VSOP image of 1156+295 at 
1.6 GHz (Fig. 3c) is comparable to that of ground based 
VLBI images at 5 GHz. This allows us to study the source 
structure on the same scales but at different frequencies. 
We note that a similar curved jet structure is visible in 
both the 1.6 GHz VSOP and 5 GHz ground-based images. 

The lower resolution VLBA image of 1156+295 at 
1.6 GHz (Fig. 3h) also displays a curved jet. In particular, 
the jet is observed to turn sharply at a few tens of milliarc- 
seconds from the core, aligning with the direction of the 
kpc jet. A high degree of linear polarization is detected in 
the 1.6 GHz VLBA data: two distinct polarized compo- 
nents are observed in the areas where the jet bends (see 
Fig. 4). A peak polarized brightness of 14.4 mJy/beam was 
detected in the core. The strongest polarized component is 
located at 2.6 mas north of the core with a peak polarized 
brightness of 18 m Jy/beam. The secondary polarized com- 
ponent was detected at 8.1 mas from the core at P.A.=:60° 
with a peak polarized brightness of 8.5 m Jy/beam. The 
percentage polarization of the core, the strongest compo- 
nent and secondary component are 1.5%, 2.3%, and 7.7%, 
respectively. The strongest polarized jet component and 
the secondary polarized jet component have perpendicu- 
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Fig. 2. A MERLIN 1.6 GHz image at the epoch 1997.41 (a) and 8.5 GHz VLA image at the epoch 2000.9 (b). 
The images were restored with a circular beam of 250 mas (FWHM) in diameter. Peak flux densities are 1.55 and 
1.0 Jy/bcam, the r.m.s. noise levels are 0.4 and 0.13 mJy/beam and the lowest contours are 1.5 and 0.4 mJy/beam, 
respectively. Contour levels increase by a factor of 2. 



Table 3. The parameters of the images in Fig. 3 



No. 


Epoch 


Array 


Band 


a* 

^peak 


r.m.s. 




Contours 




Restoring 


Beam 








(GHz) 


(Jy/beam) 


(mJy/beam ) 




(mJy/beam) 




MajxMin(mas 


), P.A 


a 


1996.43 


VLBA 


5.0 


2.08 


0.5 


1.5x( 


-1,1,2,4,..., 


1024) 


3.7 X 1.2, 


-1.8 


b 


1997.14 


EVN 


5.0 


1.02 


0.6 


2.0 x( 


-1,1,2,4,..., 


256) 


3.0 X 2.6, 


77 


c 


1997.42 


VSOP 


1.6 


0.89 


1.0 


3.0x( 


-1,1,2,4,..., 


256) 


4.4 X 1.15, 


-5 


d 


1998.12 


Global 


5.0 


0.38 


0.13 


0.5x( 


-1,1,2,4,..., 


512) 


2.3 X 0.78, 


-7.7 


e 


1999.01 


VLBA 


15.0 


2.22 


0.1 


0.3x( 


-1, 1,2,4, 


4096) 


1.1 X 0.63, 


-1.7 


f 


1999.14 


EVN 


5.0 


1.34 


0.25 


1.0x( 


-1,1,2,4,..., 


1024) 


2.3 X 2.0, 


4.2 


g 


1999.45 


Global 


5.0 


1.74 


0.22 


0.7x( 


-1,1,2,4,..., 


2048) 


2.5 X 0.85, 


-1.7 


h 


2000.15 


VLBA 


1.6 


1.0 


0.12 


0.5x( 


-1,1,2,4,..., 


1024) 


8.3 X 5.2, 


-5.5 



* peak brightness 



lar E- vectors to each other (Fig. 4). Highly polarized com- 
ponents appear to be associated with the sharp bends of 
the jet. This could be caused by a change of opacity along 
the line of sight or a transition from the optically thick to 
the optically thin regime. 

The VLBI images of the source 1156+295 reveal a core- 
jet structure with an oscillatory jet on mas scales. Two 
sharp bends in the opposite directions occur (one curves 
anti-clockwise to the north-east within a few mas of the 
core, while the other bends clockwise to the north-west 
within a few tens of mas from the core). The VLBI jet 
then aligned with the direction of the MERLIN jet. The 



oscillations in the jet structure on mas scales resemble a 
3D helical pattern. 

3.3. Light Curves at Radio Frequencies 

Fig. 5 shows the hght curves of 1156+295 at 4.8, 8.0, and 
15 GHz measured at the University of Michigan Radio 
Astronomy Observatory (UMRAO). Several flares from 
1980 to 2002 were detected. 

The intensity of a flare is stronger at higher fre- 
quencies. Flares at different frequencies do not oc- 
cur simultaneously: high frequency flares appear early. 
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Fig. 3. The VLBI Images of 1156+295. The parameters of the images are Usted in Table|31 



Similar results have been reported in some other sources The strongest flare recorded over the last 20 years 
l|Bower et al. 19971 IWehrle et al. 19981 IZhou et al. 2n00| . started at the beginning of 2001. The available data show 
The lags are consistent with the fact the flares components a maximum in the flux density of 5.7 Jy at 14.5 GHz 
are optically thick and short wavelengths peak first, which in September 2001. The next strongest flare began in 
can be seen directly from the light curves (Fig. 5). October 1997 at 8.4 and 14.5 GHz, and reached its peak 

in October 1998 with a flux density of 3.6 Jy at 14, .5 GHz. 
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Fig. 5. Light curves of 1156+295 at 4.8, 8.0 and 15 GHz from the monitoring data of the University of Michigan Radio 
Astronomy Observatory. 





Fig. 4. The VLBA image at 1.6 GHz with 

total intensity contours at 0.5 mJy/beam Fig. 6. A helical pattern for 1156+295: vector S is the 

X (-1, 1, 2, 4, 8, 16, 64, 256, 1024), and superposed sticks direction of the line of the sight, vector V is the direction 

show the orientation of electric vectors (polarization line of the velocity of the jet, 0{t) is the viewing angle between 

1" = 1.92 Jy/beam). V and S. 



Another similar peak occurred in September 1999. Such 
double peaked flares are also seen in the period from 1985 
to 1987. The separations of the double peaks are about 
one year. 

A possible explanation for the 'double-peak' is that 
only one jet component is ejected with a helical trajectory 
on the surface of a cone. The emission of the jet compo- 



nent is then enhanced along the path but only when it 
points toward the observer. If two nearby pitch angles (or 
multi pitch angles) are close to the line of the sight, we can 
observe double peaked flares (or multi peaked flares) with 
the maximum Doppler boosting from one ejected compo- 
nent. In the case of 1156+295, this results in two distinct 
peaks in the total flux density measurements. If the jet 
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Table 4. The results of model fitting at kiloparsec scales 



Observation inf. 


Stotal 


Comp 


Sc 


r 


P.A. 


a 


b/a 


P.A. 




(mjy) 




(mJy) 


(mas) 


(°) 


(mas) 




(°) 


(1) 


(2) 


(3) 


( 4) 


(5) 


(6) 


(7) 


(8) 


(9) 


1997.14 


1668 


D 


1608 


0.0 


0.0 


8.4 


0.25 


15.2 


MERLIN 




D6 


21.5 


53.7 


-8.5 


18.6 


1.0 





5 GHz 




D4 


11.4 


515 


-19.3 


132 


1.0 









D3 


13.2 


737 


-18.2 


76.9 


1.0 









D2 


7.4 


1027 


-10.8 


101 


1.0 









Dl 


13.9 


1948 


-18.4 


135 


1.0 





Errors 






8% 


9% 


3 


9% 




3 


1999.14 


2290 


D 


2214 


0.0 


0.0 


8.6 


0.38 


15.6 


MERLIN 




D6 


28.5 


56.9 


-6.0 


23.9 


1.0 





5 GHz 




D4 


13.1 


545 


-20.6 


130 


1.0 









D3 


14.0 


743 


-18.0 


70.0 


1.0 









D2 


11.2 


1057 


-12.0 


138 


1.0 









Dl 


16.6 


1938 


-18.3 


128 


1.0 





Errors 






8% 


9% 


3 


9% 




3 


2000.92 


1666 


D 


1640 


0.0 


0.0 


5.8 


0.6 


25.1 


VLA 




D3 


7.5 


747 


-16.1 


102 


1 





22.5 GHz 




Dl 


4.6 


1960 


-17.9 


88.0 


1 





Errors 






7% 


7% 


3 


7% 




3 


1997.41 


1768 


D 


1540 


0.0 


0.0 


10.1 


0.7 


15.0 


MERLIN 




D5 


27.8 


201 


-27.6 


46.9 


1.0 





1.6 GHz 




D4 


15.0 


480 


-19.7 


70.1 


1.0 









D3 


42.1 


716 


-18.4 


105 


1.0 









D2 


17.4 


1109 


-12.8 


183 


1.0 









Dl 


37.6 


1945 


-18.7 


137 


1.0 









DN 


45.7 


2055 


-2.0 


1409 


0.5 


59.5 






Db 


20.6 


138 


147.6 


156 


1.0 









Da 


11.1 


1398 


144.9 


694 


1.0 









DS 


27.6 


2152 


-161.5 


1217 


0.6 


-75 


Errors 






8% 


9% 


3 


9% 




3 


2000.92 


1042 


D 


1000 


0.0 


0.0 


13.5 


0.7 


22.3 


VLA 




D5 


2.2 


198 


-26.5 


70.4 


1.0 





8.5 GHz 




D4 


3.2 


507 


-22.9 


85.0 


1.0 









D3 


10.1 


734 


-17.6 


135 


1.0 









D2 


3.8 


1128 


-13.6 


231 


1.0 









Dl 


9.2 


1945 


-18.2 


135 


1.0 









DN 


6.6 


2121 


-2.3 


1305 


1.0 





Errors 






7% 


7% 


3 


7% 




3 



components always travel along a helical path aligned with 
the line-of-sight then double peaks should always appear 
in the light curves. In particular, double peaks should ap- 
pear at the outburst near 2002. In this shown in 
Fig. 6, we find one maximum for Doppler boosting at the 
each helical cycle. Then the peak separation time of dou- 
ble peaks Atpfe should be the period of the first detectable 
helical orbit and can be directly inferred as the period of 
the precession of the jet-base. The pure geometrical model 
described here is similar to the lighthouse model proposed 
by Camenzind & Krockenbergcr (1992|. However, in the 
case of 1156-1-295, the observed period of the peak of its 
flux density is somewhat longer than that in the lighthouse 
model. 

Two jet components ejected one after the other can 
also explain the 'double-peak' pattern. If these two jet 
components move out from the radio core, they might 



be resolved from each other by high resolution VLBI im- 
ages. Further multi-epoch VLBA observations at higher 
frequencies (2cm, 7 and 3 mm) are currently underway, 
which will help to clarify the issue. 

The peak brightness detected in our multi-epoch VLBI 
observations varies in phase with the total flux density 
measurement made at UMRAO. This indicates the most 
outbursts take place in the core area and/or in the unre- 
solved inner-jet. 

4. Morphology and structural variability on 
kpc-scales 

4.1. Model fitting the MERLIN and VLA data 

The MERLIN and VLA visibility data were fitted with el- 
liptical Gaussian components using DIFMAP. The results 
of model-fitting are presented in Table 0] Column 1 gives 
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Fig. 7. The differential image of the two 5 GHz MERLIN 
data-sets from the epoch 1999.14 and 1997.14, Peak flux 
density is 0.6 Jy/beam. Contour levels are 0.5 mJy/beam 
X (-4, -2, -1, 1, 2, 4, 8, 16, 32, ...), FWHM: 80 x 80 mas. 

the observation epoch, array, and observed frequency. The 
total CLEANed flux density of the image is listed in col- 
umn 2, followed by the component's name in column 3. 
Column 4 shows the flux density of the component. This 
is followed in columns 5 and 6 by the radial distance and 
position angle of the component (relative to the core com- 
ponent). The next three entries (columns 7 to 9) are the 
major axis, axis ratio and position angle of the major axis 
of the fitted component. The uncertainties of the model 
fitting listed in Tabic ^ are estimated using the formulae 
given by Fomalont |1999|). 

Besides the core component D, six components (Dl 
to D6) are fitted in the northern jet (see Fig. lb). Two 
counter jet components (Da and Db) and two lobe compo- 
nents (DN and DS) are detected in the 1.6 GHz MERLIN 
data (as labelled in Fig. 2a). D5 is too faint to be well 
fitted in the 5 GHz MERLIN data (Fig. la&b). D6 is not 
fitted with the 1.6 GHz MERLIN and 8.5 GHz VLA data 
because of the limitation of the resolution (Fig. 2). Only 
three main components (D, Dl and D3) are detected in 
the 22.5 GHz VLA data (Fig. Ic). 

From columns 2 and 4 of Table 01 it is clear that more 
than 95 % of the flux density came from the core at all 
frequencies, except 87% at 1.6 GHz. This is consistent with 
the flux density of the core being enhanced by beaming 
effects while the extended emission is not. 

The two main kpc-scale components (Dl and D3) are 
almost at the same position angle, ^ —18°, in all MERLIN 
and VLA images, and their FWHMs are comparable at 
1.6, 8.5, and 22 GHz, while D3 is more compact than Dl 



at 5 GHz. We note that D2 and D4 sometimes appear to 
have larger FWHMs than that of Dl and D3, which may 
be attributed to using circular Gaussian components for 
model fitting. 

A counter jet component (Db) is fitted in the 1.6 GHz 
MERLIN data (see Fig. 2a). It is located at about 138 mas 
from the core at p.a. ~ 147°. Two other low brightness 
components are also detected (labelled as Da and DS in 
Table Hand Fig. 2a). 

The lobe to counter-lobe (the components DN and NS) 
flux density ratio at a distance of 2 arcseconds from the 
core is Jiobe = 1-66. The distance and size of the counter- 
jet component Da allow us to assume that its correspond- 
ing structural pattern in the jet should include the com- 
ponents D2 and D3. Under this assumption, the flux den- 
sity ratio (D2-I-D3) to counter-jet Da is Jkpc-jet = 5.36. 
No clear jet component corresponding to Db is evident in 
Fig. 2a. This might be due to the corresponding jet com- 
ponent being embedded in the core. 

4.2. Structural variability on kpc-scales 

The two 5 GHz MERLIN observations with a time sepa- 
ration of two years (the epochs 1997.14 and 1999.14) allow 
us to study structural variability. As we mentioned in sec- 
tion 3.3, the source was in its lowest state around 1998, 
immediately before a strong flare occurred (see Fig. 5). 
The model fitting results show that the flux density of the 
core component D increased by 37% from 1997 to 1999 
(Table 01 , which is most easily explained in terms of an 
outburst occurring in the inner core. 

To study the variation in the jet radio structure di- 
rectly from the images, we subtracted the image in Fig. la 
from that in Fig. lb [Flh — Fla) to obtain a differential 
image (Fig. 7). In order to minimize the adverse effects of 
a limited dynamic range in the presence of the bright core, 
we first subtracted the core component of 1.5 Jy from the 
two MERLIN data-sets (epochs 1997.14 and 1999.14), re- 
imaged the residual uv-datasets, and then produced the 
differential image presented in Fig. 7. 

The largest difference in the fiux densities comes from 
the core. The peak brightness of the differential image is 
0.6 Jy/Beam. Some variations at a level of 1 mJy per beam 
in the jet were found. 

The increased flux density in the core can be explained 
as a flare in the compact unresolved component, as we can 
see in the light curve in Fig. 5. Both positive and negative 
variations are found above the uncertain set by the noise in 
the original images (see Fig. la&b) in the first 0.5" of the 
jet (see Fig. 7). This may be attributed to the continuous 
movement of the jet. The variations in the areas of Dl and 
D3 can not be affected by the fiare from the core, since 
it would take a few thousand years to propagate through 
the distance of 2" at the speed of light. This leads to the 
conclusion that some flares may occur, independently from 
the core, in the knot (D3) and hotspot (Dl). 
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Fig. 8. (a): A simultaneous spectral index distribution between 8.5 and 22.5 GHz at the epoch 2000.92 (color), 
superimposed on the 8.5 GHz VLA image; the wedge at the bottom shows the spectral index range —2.5 to 0.5; (b): 
A 1.6 - 5 GHz non-simultaneous spectral index distribution (color), superimposed on the 5 GHz MERLIN image at 
the epoch 1997.14; the wedge at the bottom shows the spectral index range —2.4 to 0.08. The contours levels in both 
b&c are 0.5 mJy/beam x (-1, 1, 4, 16, 64, 256) with FWHM of 200x200 mas. 
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Fig. 9. Slice plots of spectral index at p. a. -18°, corresponding to Fig. 8a and Fig. 8b respectively. 

4.3. Spectral index distributions on kpc-scales 22.5 GHz is presented in Fig. 8a by aligning the peak 

emission of the cores. An upper limit to the opacity shift 
Multi-frequency observations with the VLA and MERLIN ^3 estimated for the optically thin compo- 

allow us to produce a spectral index distribution of the ^^^^ 133^ ^^lile a frequency-dependent position difference 
source on the kpc scale (S^ cx j/"). A spectral index distri- 
bution obtained with simultaneous VLA data at 8.5 and 
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Bands 


D 


Dl 


D2 


D3 


D4 


8.5 


- 22.5 GHz 


0.51 


-0.7 




-0.3 




1.6 


- 5.0 GHz 


0.04 


-0.9 


-0.78 


-1.1 


-0.25 



of the core of 7 mas was estimated with Lobanov's model 
HT.oha,nov 199S|I. 

We also present a non-simultaneous spectral index dis- 
tributions based on the MERLIN data between 1.6 GHz 
(epoch 1997.41) and 5 GHz (epoch 1997.14) in Fig. 8b 
without consideration of the opacity shift. An upper limit 
to the opacity shift of 20 mas was estimated by compar- 
ing the modelfit results of jet component D3, and 7.4 mas 
of frequency-dependent position difference of the core was 
estimated l|Lobanov 1998|l . 

The spectral index distributions of the source (Fig. 
Sa&b) show that the strongest compact component has a 
flat spectrum. 

A steep spectrum ring appears around the core in 
Fig. 8a since the size of the core at 22 GHz is only about 
half of that at 8.5 GHz. 

The kpc scale jets clearly have a steeper spectrum than 
the core (see Fig. 8a&b). There is also a general steepening 
of the radio emission towards the edges of the jet emission. 
Some isolated flat spectrum spots in the jet are also ob- 
served, especially in Fig. 8b. These may be real or due to 
the observations being made at different epochs. 

Slice plots of the spectral index distributions along the 
jet axis at position angle of -18° are shown in Fig. 9, cor- 
responding Fig. 8a and 8b, respectively. The slice lines 
do not intersect the peak of the each component since the 
components are not aligned on a straight line. We can also 
estimate the spectral indexes with the modelfit results of 
Table 3 (see Table E)). 

It is clear that the core has a flat spectral index, while 
the jet components Dl, D2, , D3 and D4 have steep spec- 
tral indexes. Each component has a steep spectrum edge. 

5. Morphology and structural variability on 
pc-scales 

5.1. Model fitting of the VLBI data 

The VLBI self-calibrated data-sets associated with each 
image in Fig. 3 were fitted with Gaussian components us- 
ing DIFMAP. An elliptical Gaussian component was fitted 
to the strongest component, and circular Gaussian bright- 
ness distributions were used for the jet components to es- 
timate their sizes. 

The parameters of model fitting are listed in Table |H1 
Column 1 gives the observation epoch, array, and fre- 
quency, columns 2 to 8 are the same as the columns 3 
to 9 in Table 0] The uncertainties of the model fitting 
are estimated in the same way as for the acrsecond-scale 
images (see section 4.1) and are listed in the table. 



I 




Right Ascension (mas) 



Fig. 10. The model image at the epoch 1999.45 at 5 GHz. 
Peak brightness is 1.71 Jy/beam. Contour levels are 
1.5 mJy/beam x (-1, 1, 1, 4, 16, 64, 256). 



Five components are fitted to the VLBI data. In order 
to compare our results with previous publications, we label 
the components following the convention introduced by 
Finer and Kinghani (1997) from C to C4 (see Fig. 10). 

Fig. 10 is the model-fitted image (Gaussian compo- 
nents restored with a 2.4 x 1.1 mas, at p. a. — 2.1° beam) 
of 1156-1-295 based on the Global VLBI data at the epoch 
1999.45 at 5 GHz. This is in good agreement with the 
CLEANed image shown in Fig. 3g. 

The core contributes a higher fraction of the total flux 
density at higher frequencies. The fraction of the core flux 
density is about 97% at 15 GHz, 80 - 86% at 5 GHz (except 
the lowest level of 65% at the minimum of activity at 
the epoch 1998.12), and 65 - 70% at 1.6 GHz. This is a 
clear indication that the high frequency emission comes 
from the inner regions of the source. The fraction of the 
core flux density is higher during active outburst phases 
(e.g. 85% at the epoch 1999.45 at 5 GHz) compared to 
more quiescent epochs (e.g. 65% at the epoch 1998.12 at 
5 GHz). This suggests that most of the flare occurred in 
an area less than 1 mas in size. 

For the individual components, we note: a) we use one 
elliptic Gaussian component CI to fit the two separated 
outer emission regions in the north-east area at epoch 
1996.43 (Fig. 3a); b) CI is not detected at the epoch 
1997.14 due to sensitivity limitations and is resolved at 
15 GHz (epoch 1999.01); c) components C2 and C3 are 
not resolved in the 1.6 GHz VLB A observations and are 
fitted as a combined centroid component C2&3; d) C4 is 
detected after 1999 at 15 GHz; e) outer component CO is 
detected at a distance of 30 mas from the core (Fig. 3h). 
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Table 6. Model-fitting of the VLBI data 



Obs. inf. 


Comp 


Sc 


r 


P.A. 


a 


b/a 


p.a. 


Seq 


Sic 


7 


9 


n 


Tr 






(mjy) 


(mas) 


n 


(mas) 




n 








n 


K 


K 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


1996.43 


C 


2097 


0.0 


0.0 


0.55 


0.24 


-30.6 


15.8 


15.5 






2.5x10'^ 


1.6x10'' 


VLBA 


C3 


141.2 


3.22 


-2.9 


1.48 


1.0 









11.4 


3.4 


5.3x10^ 




C Band 


C2 


133.9 


4.92 


18.0 


1.88 


1.0 









13.8 


3.6 


3.1x10^ 






CI 


129.5 


11.01 


34.9 


9.05 


0.67 


-14.9 










2.0x10^ 




Errors 




12% 


15% 


7 


16% 




7 














1997.14 


C 


1109 


0.0 


0.0 


0.51 


0.41 


-25 


5.0 


6.0 






0.9x10''' 


1.8x10" 


EVN 


C3 


117.8 


3.95 


-3.6 


1.69 


1.0 









13.6 


8.9 


3.5x10^ 




C Band 


C2 


172.4 


5.06 


20.2 


3.93 


1.0 









21.2 


7.5 


9.4x10® 




Errors 




18% 


20% 


10 


20% 




7 














1997.42 


C 


927.3 


0.0 


0.0 


0.90 


1.0 





5.3 


4.3 






0.9x10'^ 


1.7x10" 


VSOP 


C3 


135.0 


3.01 


7.2 


2.05 


1.0 









13.1 


8.8 


2.7x10^° 




L Band 


C2 


240.2 


5.66 


18.2 


4.11 


1.0 









20.3 


7.3 


1.2x10^° 






CI 


104.4 


10.97 


35.5 


5.73 


1.0 













2.7x10^ 




Errors 




n% 


15% 


10 


20% 




7 














1998.12 


C 


504.0 


0.0 


0.0 


1.0 


0.39 


1.4 


0.5 


0.9 






l.lxlO''' 




Global 


C3 


89.9 


3.71 


-1.6 


1.75 


1.0 













2.5x10^ 




C Band 


C2 


95.4 


5.59 


21.1 


2.66 


1.0 













1.1x10^ 






CI 


82.6 


10.10 


33.1 


7.84 


1.0 













1.1x10* 




Errors 




18% 


10% 


10 


10% 
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1999.01 


C 


2245 


0.0 


0.0 


0.15 


0.37 


15.4 


17.1 


22.4 






2.5x10^'' 


1.4x10" 


VLBA 


C4 


26.3 


1.25 


-4.6 


0.8 


1.0 









12.7 


3.3 


3.8x10® 




U Band 


C3 


17.7 


4.19 


-2.6 


1.5 


1.0 









11.8 


3.0 


7.4x10^^ 






C2 


28.0 


6.23 


16.0 


2.6 


1.0 









14.0 


3.1 


3.8x10^^ 




Errors 




15% 


13% 


5 


13% 




7 














1999.14 


C 


1387 


0.0 


0.0 


0.39 


0.40 


3.4 


12.3 


11.8 






2.0x10'^ 


1.6x10" 


EVN 


C4 


42.2 


1.55 


-9.6 


1.1 


1.0 









11.9 


4.6 


2.9x10^^ 




C Band 


C3 


31.9 


4.26 


-5.7 


1.2 


1.0 









10.6 


4.6 


l.SxlO'' 






C2 


76.3 


5.83 


15.8 


2.3 


1.0 









13.7 


4.6 


1.2x10^^ 






CI 


71.7 


10.62 


30.0 


6.6 


1.0 













1.3x10® 




Errors 




15% 


15% 


6 


15% 




7 














1999.45 


C 


1765 


0.0 


0.0 


0.35 


0.44 


-6.9 


18.5 


12.4 






2.7x10'^ 


1.4x10" 


Global 


C4 


52.9 


1.55 


-6.5 


1.33 


1.0 









13.1 


2.8 


2.4x10^ 




C Band 


C3 


54.1 


4.31 


-3.7 


2.06 


1.0 









12.3 


2.7 


l.lxlO'' 






C2 


108.5 


6.14 


20.0 


3.55 


1.0 









14.3 


3.0 


7.4x10® 






CI 


92.1 


11.2 


33.3 


8.07 


1.0 













1.2x10® 




Errors 




9% 


10% 


5 


10% 




5 














2000.12 


C 


931.3 


0.0 


0.0 


1.6 


0.28 


-19.1 


6.1 


4.8 






l.lxlO^'' 


1.8x10" 


VLBA 


C2&3 


245.1 


5.23 


16.4 


3.9 


1.0 









15.1 


7.5 


1.3x10" 




L Band 


CI 
CO 


195.0 
88.1 


10.53 
31.66 


31.2 
-1.1 


7.9 

25.4 


1.0 
1.0 














2.5x10^ 




Errors 




10% 


12% 


5 


13% 




5 















5.2. The differential VLBI image at 5 GHz 

As noted in section 5.1, the largest structural variations 
during the flare appeared in the vicinity of the core. The 
two high sensitivity global VLBI observations at 5 GHz 
(epochs 1998.12 and 1999.45) were made in the states of 
low and high flaring activity respectively. The differential 
image of these two global VLBI datasets {F3g — F3d) is 
shown in Fig. 11. Most of the different flux density in the 
image comes from the core area. The peak (residual) flux 
density is 1.3 Jy/beam, which is about 3 times of that 
of the total intensity peak observed at the epoch 1998.12. 



This large difference in the core region is clearly associated 
with the violent flare in the core that has occurred after 
the epoch 1998.12 and during the epoch 1999.45. 

The existence of a 'negative emission' component in 
the differential image (Fig. 11) can be explained by a vari- 
ation in the direction of the jet. When the direction of the 
relativistic jet comes close to the line of sight, its flux den- 
sity is enhanced by a Doppler boosting; when the orienta- 
tion of jet changes from the line of sight, the component 
becomes fainter. This pattern is characteristic of a rela- 
tivistic jet with a helical trajectory. 
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Fig. 11. The differential image of the two 5 GHz Global 
VLBI data-sets from the epochs 1999.45 and 1998.12, 
the peak brightness is 1.3 Jy/beam., contour levels are 
2.0 mJy/beam x (-8, -4, -2, -1, 1, 2, 4, 8, 16, 32, 512). 
FWHM: 2 x 1 mas, at 0° 




5.3. Proper motion of the jet components 

The positions of the components during the period from 
1996 to 2000 are shown in Fig. 12. The x-axis is the time 
and the y-axis shows the distance of the components from 
the core. Different symbols are used to represent different 
observing frequencies: a diamond for 1.6 GHz, a star for 
5 GHz, and a square for 15 GHz. We also plot a triangle 
for the 22 GHz component labelled as B2 by Jorstad et al. 
(2001). We identify B2 is the same with the component C4 
in this paper (see Table 4, Fig. 10). The components C2 
and C3 are not resolved at 1.6 GHz at the epoch 2000.12. 
This combined centroid is shown as a filled diamond sym- 
bol. 

Since the positions of the components could be fre- 
quency dependent, we estimate the proper motion of com- 
ponents C2, C3 and C4 via the 5 GHz data sets. These 
are 0.39±0.1, 0.30 ± 0.08 and 0.34 ± 0.08 mas/yr respec- 
tively. These values correspond to the apparent velocities 
of 13.7 ± 3.5, 10.6 ± 2.8, and 11.8 ± 2.8c. 

The component CI is located at the point where the jet 
turns sharply. It does not appear to have any appreciable 
motion in the r direction, while the position angle changed 
from about 35°to 31°, which correspond to a apparent 
velocity of about 7c. 

Jorstad et al. (2001) reported the proper motion mea- 
sured at 22 GHz of three components in this source, the 
period of the VLBA observation is overlapped with ours, 
the reported apparent velocity of B2 (labelled as C4 in 
this paper) is 11.8 ± 1.2 (corrected to go — 0.5), based 
on three epochs during period 1996.60 to 1997.58, which 
is in agreement with our result. The B3 component in 
their paper is unresolved in our VLBI data due to the 
limited resolution at low frequencies. Finer & Kingham 
(1997) reported that the apparent velocities of CI to C4 
are 13.5 ± 3.5, 8.1 ± 1.7, 8.5 ± 1.4, and 5.4 ± 1.8 c, respec- 
tively, based on their geodetic data obtained during the 
period 1988.98 to 1996.25. 

The apparent radio component velocity of 40c is the 
highest super-luminal velocity reported for AGN to date 
HMcHardy et al. 1990||McHardy et al. 1993| ). However, as 
this value was derived based on four-epoch of VLBI ob- 
servations at three frequencies, the apparent motion could 
be 'caused' by frequency-dependent effects. 



Fig. 12. Apparent proper motion of the jet components: 
1.6 GHz (diamond), 5 GHz (star), 15 GHz (square), 
22 GHz (triangle, Jorstad, et al. 2001). 



The observed feature can be also explained directly 
by the proper motion of the jet. If one component moves 
from point A to point B, an observer will see a negative 
component at the point A and a positive component at 
point B with the same absolute value of flux density. This 
seems not like in the case of Fig. 11. 



5.4. Spectral index distribution on pc-scales 

Although no simultaneous VLBI data of the source at dif- 
ferent frequencies are yet available, we attempted to esti- 
mate the spectral index distribution of the jet on pc scales 
using the multi-epoch/multi-frequency images presented 
here. The non-simultaneous spectral index distribution be- 
tween 1.6 (epoch 2000.15) and 5 GHz (epoch 1999.45), as 
wefl as between 5 (epoch 1999.45) and 15 GHz (1999.01) 
are presented in Fig. 13 by summing the zero opacity shift 
in the cores. The upper limit opacity shift of 0.3 and 0.7 
mas were estimated by comparing the modelfit results of 
the optical thin jet component for Figs. 13a and 13b, re- 
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Fig. 13. (a): The 5 vs 15 GHz spectral inJ^S^'distribution (color), superimposed on the 15 GHz VLB'A'Image at the 
epoch 1999.01, the wedge at the bottom shows the spectral index range —1.95 to 0.25; (b): the 1.6 vs 5 GHz spectral 
index distribution (color), superimposed on the 5 GHz VLB A image, the wedge on the bottom shows the spectral 
index range —1.5 to 0.62. The contours levels in a & b are 2 mJy/beam x(— 1, 1,4, 16,64,256) with FWHM of 1.5x1 
mas and 7x4 mas, respectively. 



Table 7. The spectral index on the pc scale 



Bands 


C 


CI 


C2 


C3 


C4 


5 - 15 GHz 


0.22 




-1.2 


-1.0 


-0.6 


1.6 - 5 GHz 


0.58 


-0.7 









spectively. Meanwhile, 0.5 mas of the uncertainty of the 
frequency-dependent position difference of the core was 
estimated with Lobanov's model IjLobanov 1998|l for the 
spectral index maps Figs. 13a and 13b. 

Fig. 13a shows the 5—15 GHz spectral index distri- 
bution. It is clear that the core area has a flat spec- 
trum a ~ 0.25. The emission becomes a steep spectrum 
{a ^ —0.7 to —1.5) outwards from the core. 

Fig. 13b is the 1.6—5 GHz spectral index distribution. 
There is an inverse spectrum (a ~ 0.58) in the core area. 
Two flat spectrum features are seen in the jet. 

Figs. 14a and 14b show the profile of the 5—15 GHz 
spectral index along the direction of P.A.=0° and 1.6 — 
5 GHz spectral index along P.A.= 25°, corresponding to 
Figs. 13a and 13b, respectively. 

Table [71 gives the spectral index values estimated from 
the model-fit components listed in Table [7| 

6. Physical parameters of the core and jet 

6.1. The Equipartition Doppler factor 5eq 

As we mentioned in the previous sections, the flares of 
1156-1-295 may be enhanced by the Doppler boosting. 



The Doppler factor of the outflow from compact radio 
cores can be estimated assuming the energy equiparti- 
tion between the radiating particles and the magnetic field 
IjReadhead 19941 |Uuij osa fc Daly 1996| ). 

The method was considered using the radio observed 
parameters at the turnover frequency, but actual VLBI 
observations were not done at the turnover frequency. We 
will use the radio parameters we obtained from our VLBI 
observations to roughly estimate the Doppler boosting in 
1156-^295. 

By following the assumptions of Scott & Readhead 
(1977) and Marscher (1987), the Equipartition Doppler 
factor 6eq of the core component for each epoch is esti- 
mated and listed in column 9 of Table |H| The values are 
in the range of 5 to 18.5 except 0.5 at the epoch 1998.12. 
The latter is expected, since the source is at its lowest 
state. 

6.2. Inverse Compton Doppler factors Sjc 

For comparison, the Doppler factor of the core component 
can also be estimated with the Synchrotron self-Compton 
emission in a uniform spherical model l)Marscher 1987|l . 

Under the assumption in section 6.1, we estimated Sjc 
using the equation (1) of Ghisellini et al. (1993) with 
the 0.075 /iJy X -ray flux density of 1156+295 at 2 keV 
dMcHardy 1985) ). The inverse Compton Doppler factors 
Sic of 1156+295 at all epochs were listed in the column 
10 of Table^l The values cover a wider range from 4.3 to 
22.4 except 0.9 at the epoch 1998.12. 
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Fig. 14. (a): Slice plot of the spectral index at p. a. 0° corresponding to Fig. 13a; (b): Slice plot of the spectral index 
at p. a. 25° corresponding to Fig. 13b 



The estimated values may be a lower limit since part of 
X-ray emission could come from the thermal or other emis- 
sion. The variability of X-ray flux and non-simultaneous 
observations between radio and X-ray are causes of uncer- 
tainty in estimating Sjc- 

The Sic and the S^q are comparable. The average val- 
ues of Sjc is about 96% of that of Seq- The Doppler factor 
varied with time. The variation was correlated with the 
total flux density variation of the source. The core has 
higher Doppler factor at a higher frequency than that at 
lower frequencies. 

6.3. Lorentz factor and viewing angle 

In the relativistic beaming model (Sapp depends on the 
true f3{v/c) factor and the angle to the line of sight 9 
) |Rees fc Simon, 1968, ), 



Papp — r) n ■ 

1 — P COS 

The Doppler factor can be written as 

(5 = 7^1(1 -/3cos6i)-^ , 



(1) 



(2) 



where 7 = 1/^1 — (3'^. 

The Lorentz factor 7 and the viewing angle ip can be 
computed (Ghisellini et al. 1993): 
j;2 



7 



2 +5'' 

app ' 



1 



26 



tan (p 



2/3, 



app 



PI, 



(52 _ 1 



(3) 

(4) 



The 7 and (f) of C2, C3, and C4 shown in columns 11 
and 12 of Table 4 are estimated from the Doppler factor 



{5 — Seq is used) and their apparent velocities f3app- The 
Lorentz factor of each epoch changes between 13.8 and 
21.2, 10.6 and 13.6, 11.9 and 13.1 for C2, C3 and C4, 
respectively, while the viewing angles varied in the range 
of 2.7° to 8.9°. 

Ghisellini et al. reported that the jet to counter-jet 
brightness ratio can be estimated from (1993): 



l + /3cos6> 



(5) 



Taking Jkpc~jet = 5.36 and Jiobe = 1-66. (see section 
4.1), and assuming the spectral index of the kpc jet and 
the lobe a = —0.9 (Tablc|SJ), we have: 



{P cos 9) kpc- jet = 0.3 
(/3cOS0)iobe =0.1 . 



(6) 
(7) 



It indicates that the kpc jet is faster and/or has a 
smaller viewing angle than the lobe components. 

6.4. The brightness temperature 

A high brightness temperature in flat spectrum radio 
sources is usually considered to be an argument sup- 
porting relativistic beaming effect. In the synchrotron 
model, with the brightness temperature close to ~ 10^^ K, 
the radiation energy density dominates the magnetic 
energy density, and produces a large amount of en- 
ergy in a very short time scale. The relativistic par- 
ticles quickly cool (Compton catastrophe). The maxi- 
mum brightness temperature that can be maintained is 
about ~ 10^2 K in the rest frame of the emitting plasma 
( jKellermann fc Pauliny-Toth 1969| ). 
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Fig. 15. The brightness temperature distribution of VLBI 
components: 1.6 GHz diamond, 5 GHz circle, 15 GHz tri- 
angle. 

The brightness temperature Tf, of an elliptical 
Gaussian component in the rest frame of the source is 
IShen et al. 199711 

Tb = 1.22 X 10i'^^(l + z)K , (8) 

where z is the redshift, 5*,^ is the observed peak flux density 
in Jy at the observed frequency i^ob in GHz, a and b are 
the major and minor axes in mas. 

The brightness temperatures of the the components 
of 1156+295 are calculated with equation (8) at various 
epochs. The results are listed in column 13 of Table 4. We 
present the brightness temperature distribution in Fig. 15. 
It is clear that the component brightness temperatures are 
inversely proportional to their observed frequencies. The 
brightness temperature of the core component is around 
10^^ K, while the brightness temperatures of the jet com- 
ponents are much lower (around 10^~^" K). 

By considering the Doppler boosting of the source, the 
intrinsic brightness temperature T,. can be estimated as 
Tr — Tb/S and listed in column 14 of Table 4 {6eq is used). 
The intrinsic brightness temperatures of the core compo- 
nent at all epochs are almost constant {r^ lO^^K). This 
may serve as an upper limit of the brightness temperature 
in 1156+295 in the rest-frame of the source. Since the 
estimated Seq is an upper limit, the intrinsic brightness 
temperatures here are the lower limit. 

7. A helical pattern 

The helical jet model for AGNs was proposed (e.g. 
[Conway fc Murphy, 19931 ) to explain a bi-modal distribu- 
tion of the difference in the orientation of arcsecond and 
mas structures in core-dominated radio sources. Helical 
jets are a natural result of precession of the base of the 



Fig. 16. Measured positions of all components at all fre- 
quencies and all epochs for 1156+295. 

jet (e.g. Linfield 1981) or fluid-dynamical instabilities in 
the interface between the jet material and the surround- 
ing medium (Hardee 1987). Assuming the initial angular 
momentum of the jet is due to the precession of the accre- 
tion disk, the helical structure observed might be related 
to the precession of a wobbling disk while the interface in- 
stability produces a wiggling structure that reflects a wave 
pattern fe.g. lZhao et al. 1992|l . 

As we mentioned in section 3.3, the trajectory of the 
jet in 1156+295 could be a helix on the surface of a cone 
(Fig. 6). Fig. 16 shows the measured position angles of all 
components at all frequencies and all epochs for 1156+295. 
The two components (B2 and D) from the Jorstad et al's 
paper (2001) are also included. The x-axis shows the dis- 
tance of the components from the core in logarithmic scale 
and the y-axis shows the position angle. This symmetric 
distribution of the position angle suggests that the jet in 
1156+295 is spatially curved and the jet components fol- 
low a helical path. 

By assuming the conservation laws for kinetic energy 
and momentum, the helical model parameters can be de- 
rived following (,StefFen et al 19951 : 

z{t)=l3,t + zo , (9) 
r{t) = I3^t i&nil} + 7-0 , (10) 

^W = 0o + --^ , (11) 

(3^ tan-f/) To 

where z,r,(j) are the cylindrical coordinates; /3 is the ve- 
locity of the jet component in units of the speed of light; ip 
is opening half-angle of the cone; t is the time in the rest 
frame of the source (jet); and lu is the angular rotation 
velocity {d(j)/dt). 

The estimated viewing angles 9 of inner jet compo- 
nents in 1156+295 are 3° to 8° (see table 4). A value of 
5° was adopted as the angle between the line of sight and 
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Fig. 17. The positions of jet components and the projec- 
tive path of a hehcal jet in the image plane on the pc-(/e/t) 
and the kpc-{right) scales. 

the axis of the cone; the half-opening angle of the cone 
'ip= 2.5° was estimated from the observed oscillatory am- 
plitude of position angle. Figs. 16 and 17 show one possible 
approximation rather than the best fit to the data. Fig. 
16 represents the relationship between the position angle 
and the projected distance from the core of this helical 
jet. Fig. 17 shows the projected path of this helical trajec- 
tory in the image plane on the two different scales. These 
figures show that a simple helical trajectory can well rep- 
resent the radio structure from 1 mas to 1 arcsecond (D2 
component). 

Beyond the D2 component, the helical path is proba- 
bly interrupted or it changes direction due to the inter- 
action between the jet and the medium. The model does 
not consider scales less than 1 mas because of the lack of 
observational data. Jorstad et al. (2001) pointed out that 
the trajectories of components can fall anywhere in the re- 
gion of position angle from —25° to 25° out to a distance 
of ^ 1 mas from the core based on their data at 22 GHz. 
This does not contradict the helical model presented here. 

Fig. 18 shows that the components move along curved 
path. For comparison, we also included the data from the 
geodetic VLBI observation by Finer & Kingham (1998) 
(estimates from their Fig. 10) during period from 1988 to 
1996 in Fig. 18, labelled as PKCl, PKC2 and PKC3. The 
two data sets show the movement of jet components along 
a curved trajectory on the time scale of 20 years. 

In addition, we note some small structural oscillations 
observed on the mas scale (e.g. components C3 and CI). 
It is not clear whether this is a processing on a short time 
scale or measurement errors. 

If the jet moves on the surface of a cone along a helical 
trajectory, the helical wavelength (A, the pitch distance) 
of this helical trajectory increases with the time. The ap- 
parent period of this helical jet (starting at about 1 mas 
from the core) at the projected distance is about 1200 
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Fig. 18. The evidence for the movement of the compo- 
nents along a curved path. 

mas. If we assume that the averaged proper motion (/i 
= 0.3 mas/yr) is constant on the 1 mas to 1 arcsecond 
scale (as we assumed), the lower limit of the helical pe- 
riod in time in the source-rest frame is about 10^ years 
(P = {1 + z) * X/ fi). Considering the possible deceleration 
of the jet on the large scale, the real period is likely to be 
longer. 

The helical jet may be due to a binary black hole in the 
nucleus of the source, such as in the BL Lac Object OJ 287 
(e.g. iSillanpaa et al. 19 88 Lchto & Valtoncn 199^. The 
model for both the optical and radio light curves ob- 
served in OJ 287 suggests that the helical jet and 
double peaked cyclic outbursts can be produced by a 
pair of super-massive black holes in a binary system 
IjVillata et al. 1998|l if the twin jets ejected from each of 
the BHs were in a precession around the orbital axis of 
the binary and the Doppler boosting effect is important. 

Alternatively, in a spinning super-massive black hole 
system, if the disk axis offsets from the spin axis, the in- 
trinsic torque of the spin black hole placing on an accre- 
tion disk could also produce a precession of the axis of the 
accretion disk l|Liu & M cha 2002J. The variation of the 
direction in the jet motion produced from a wobbling disk 
can also be responsible for the helical pattern. 

8. Summary 

We investigated the radio structures of the radio-loud 
quasar 1156-1-295 in detail on various angular scales using 
different interferometric arrays at different frequencies. 

1. The source shows an asymmetric double structure 
at 1.6 GHz with 250 mas resolution (about 1.5 kpc), the 
northern jet is relativistically Doppler boosted, the south- 
ern emission is resolved into several low brightness regions 
(see the color figure of the cover). 

2. At 80 mas resolution (about 0.5 kpc), the northern 
jet has an almost straight structure at P. A. = —18° with 
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some sinusoidal fluctuations. The source is dominated by 
a flat core, and the jet emission has a steep spectrum. 

3. On the VLBI scale, the jet bends from the north 
to the north-east at 3~4 mas from the core, and then 
turns about 90° to the north-west, thus aligning with the 
direction of the arcsecond-scale jet. The 1.6 GHz VLBA 
image links the VLBI emission and arcsecond emission of 
the jet, which follows an S type structure. This may be an 
evidence of a helical jet. The spectral index distribution 
at the mas scale shows a flat spectrum VLBI core and 
a steep spectrum jet. The difi'erential image demonstrates 
that the flux density variation mainly comes from the very 
compact core area. 

4. The proper motions in three components were de- 
tected; the apparent superluminal velocities are in the 
range of 10.6 to 13.7 c. The differential images also show 
evidence of bulk motion of the emitting material along a 
curved trajectory. 

5. Ejection of the VLBI components corresponds to 
the flares in the light curves. This is consistent with the 
hypothesis that total flux density variations mainly occur 
in the core region. 

6. High Doppler factors (^18), and a high Lorentz fac- 
tor (^^21) were measured in the northern jet components 
on pc-scales supporting the assumption that the jet moves 
relativistically. 

7. We found that a helical trajectory along the surface 
of a cone can well represent the radio structure from 1 mas 
to 1000 mas, based on the conservation laws for kinetic 
energy and momentum. The period of the helical path is 
longer than 10^ years. 

The study provides the most comprehensive account 
so far of kpc- and pc-scale morphological properties of the 
radio emission in the AGN 1156-1-295. These properties 
are analyzed in the framework of the jet helical model e.g. 
a binary black hole in the nucleus of the source. However, 
other models, e.g. based on the Kelvin-Helmholz insta- 
bilities (JIardcc 1987J, the trailing shock model in a rel- 
ativistic jet HGomez et al. 200T|) or MHD processes could 
be verifled in future studies against the observing data 
presented. 

Acknowledgements. This research was supported by the 
National Natural Science Foundation of PR China (19973103), 
and Chinese fund NKBRSF (G1999075403). The authors 
thank Dr. Jim Ulvestad for his work on the fringe-fitting and 
imaging the VSOP in-orbit checkout data of 1156-1-295. XYH 
thanks JIVE and ASTRON for their hospitality during his visit 
to the Netherlands for the data reduction in 1999 and 2002. 
The studies reported in this paper have been supported by 
the grant for collaborative research in radio astronomy of the 
Dutch (KNAW) and Chinese Academies of Sciences (CAS). 
XYH thanks Dr. Chris Carilli for the help of calibration for the 
polarization VLBA data and Dr. Jun-Hui Zhao for the useful 
discussion on the wobbling disk. The authors are grateful to the 
staff of the European VLBI Network radio observatories and 
the data reduction center, the MERLIN, the NRAO for support 
of the observations. The European VLBI Network is a joint 
facility of European, Chinese, South African and other radio 



astronomy institutes funded by their national research coun- 
cils. The National Radio Astronomy Observatory is a facility 
of the National Science Foundation operated under a cooper- 
ative agreement by Associated Universities, Inc. The authors 
gratefully acknowledge the VSOP Project, led by the Institute 
of Space and Astronautical Science in cooperation with the 
world-wide network of agencies, institutes and facilities. We 
thank the referee for critical reading of the manuscript and 
many valuable comments. This research has made use of data 
from the University of Michigan Radio Astronomy Observatory 
which is supported by funds from the University of Michigan. 

References 

Antonucci, R. R. & Ulvestad, J. S., 1985, ApJ, 294, 158 
Blandford, R. D. & Rees, M. J., 1974, MNRAS, 165, 395 
Bower, G. C, Backer, D. C, Wright, M., Forster, J. R., AUer, 

H. D. & AUer, M. F., 1997, ApJ, 484, 118 
Conway, J. & Murphy, D. W., 1993, ApJ, 411, 89 
Camenzind, M. & Krockenberger, M., 1992, A&A, 255, 59. 
Field, G. B. & Rogers, A. D., 1993, ApJ, 403 94 
Fomalont, E. B., 1999, in Synthesis Imaging in Radio 

Astronomy II, Taylor G. B., CariUi C.L., & Perley R.A. 

(eds.), P301 

Fomalont, E. B., Frey, S., Paragi, Z., Gurvits, L. I., Scott, W. 
K., Taylor, A. R., Edwards, P. G., & Hirabayashi, H., 2000, 
ApJS, 131, 95 

Garrington, S. T., Garrett, M. A., & Polatidis, A., 1999, New 

Astronomy Review, 43, 629. 
Garrington, S. T., Muxlow, T. & Garrett, M. A., in Galaxies 

and their constituents at the highest angular resolutions, 

lAU symp., 205, Eds. Schilizzi, R.T., Vogel, S.N., Paresce, 

F., et al., P102. 

Ghisellini, G., Padovani, P., Celotti, A., & Maraschi, L., 1993, 
ApJ, 407, 65 

Glassgold, A. E., Bergman, J. M., Huggins, P. J., Kinney, A. 

L., Pica, A. J., Pollock, T. J., et al., 1983, ApJ, 274, 101 
Gomez, J.-L., Marscher, A. P., Alberdi, A., Jorstad, S. G., and 

Agudo, I., 2001, ApJ, 561, L161 
Guijosa, A. & Daly, R. A., 1996, ApJ, 461, 600 
Gurvits, L. I. Kellermann, K. I., Fomalon, E. B., Zhang, H. Y., 

2003, in preparation 
Hardee, P. E., 1987, ApJ, 318, 78 

Hartman, R. C, Bertsch, D. L., Bloom, S. D., Chen, A. W., 
Deines- Jones, P., Esposito, J. A., et al. 1999, ApJS, 123, 
79 

Hirabayashi, H., Hirosawa, H., Kobayashi, H., Murata, Y., 
Edwards, P. G., Fomalont, E. B. et al., 1998, Science, 281, 
1825 

Hong, X. Y., Jiang, D. R., Gurvits, L. I., Garrett, M. A., 
Schilizzi, R. D., & Nan, R. D., 1999, New Astronomy 
Review, 43, 699 

Hong, X. Y., et al. 2003, in preparation. 

Jorstad, S. G., Marscher, A. P., Mattox, J. R., Wehrle, A. E., 
Bloom, S. D., & Yurchenko, A. V., 2001, ApJS, 134, 181 

Kellermann, K. I., Vermeulen, R. C, Zensus, J. A., Cohen, M. 
H., & West, A., 1999, New Astronomy Review, 43, 757 

Kellermann, K. I. & Pauliny-Toth, I. I. K., 1969, ApJ, 155 

Lehto, H. J., & Valtonen, M. J., 1996, ApJ, 460, 270 

Linfield, R. P., 1981, ApJ, 250, 464 

Liu, S. & Melia, F., 2002, ApJ Lett, 566, L77 

Lobanov, A. P. 1998, A&A, 330, 79 



X. Y. Hong et al.: A relativistic helical jet in the 7-ray AGN 1156+295 



19 



Marscher, A. P. 1987, in Superluminal Radio Sources, Zensus, 

J. A., & Pearson, T. J. (eds.), Cambridge University Press, 

New York, P208 
Mchardy, I., 1985, Space Sci. Rev. 40, 559 
McHardy, I. M., Marscher, A. P., Gear, W. K., Muxlow, T., 

Lehto, H. J., & Araham, R. G., 1990, MNRAS, 246, 305 
McHardy, I. M., Marscher, A. P., Gear, W. K., Muxlow, T., 

Lehto, H. J., & Araham, R. G., 1993, MNRAS, 261, 464 
Mukherjcc, R., Bortsch, D. L., Bloom, S. D., Dingus, B. L., 

Esposito, J. A., Fitchcl, R.C., ot al, 1997, ApJ, 490, 116 
Piner, B. G. & Kingham, K., 1997, ApJ, 485, L61 
Piner, B. G. & Kingham, K., 1998, ApJ, 507, 706 
Readhead, A. C. S., 1994, ApJ, 426, 51 
Rees, M. J. & Simon, M., 1968, ApJ, 152, 145 
Roos, N., Kaastra, J. S. & Hummel, C. A., 1993, ApJ, 409, 130 
Scott, M. A., & Readhead, A. C. S., 1977, MNRAS, 180, 539 
Shen Z.-Q., Wan, T.-S., Moran, J. M., Jauncey, D. L., 

Reynolds, J. E., Tzioumis, A. K., et al. 1997, AJ, 114, 1999 
Shepherd, M. C., Pearson, T. J., & Taylor, G.B., 1994, BAAS, 

26, 987 

Sillanpaa, A., Haarala, S., Valtonen, M. J., et al. 1998, ApJ, 
325, 628 

Steffen, W., Zensus, J. A., Krichbaum, T. P., Witzel, A., & 

Qian, S. J., 1995, A&A, 302, 335 
Thomasson, P., 1986, QJRAS, 27, 413 

Thomposon, D. J., Bcrtsch, D. L., Dingus, B. L., Esposito, J. 

A., Ettennc, A., Fichtcl, C. E., ct al, 1995, ApJS, 101, 259. 
Tornikoski, M. & Lahtccnmaki, A., Proceedings of the fifth 

Compton Symposium, 2000. Edited by Mark L. McConnell 

and James M. Ryan AIP Conference Proceedings, Vol. 510., 

p.377. 

Veron-Cetty, M.-P. & Veron, P., 1998, A Catalogue of Quasars 
and Active Galactic Nuclei (8th edition), ESO Sci. Report 
No. 18. 

Villata, M., Raiteri, C. M., Sillanpaa, S., & Takalo, L. O., 1998, 

MNRAS, 293, L13 
Wehrle, A. E., Plan, E., Urry, C. M., Maraschi, L., McHardy, 

I. M., & Lawson, A. J., 1998, ApJ, 497,178 
Wills, B. J., Wills, D., Breger, M., Antonucci, R. R., & 

Barvainis, R. E., 1983, ApJ, 274, 62 
Wills, B. J., Pollock, J. T., AUer, H. D., AUer, M. F., Balonek, 

T. J., Barvainis, R. E., et al., 1992, ApJ, 398, 454 
Zhao, J.-H., Burns, J. O., Hardee, P. E., & Norman, M. L., 

1992, ApJ, 387, 83 
Zhou, J. F., Hong, X. Y., Jiang, D. R., & Venturi, T., 2000, 

ApJ Letter, 540, L13 



